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A form of an indirect Ruderman-Kittel-Kasuya-Yosida (RKKY)-like coupling between magnetic 
on-site impurities in armchair graphene nanoribbons is studied theoretically. The calculations are 
based on a tight-binding model for a finite nanoribbon system with periodic boundary conditions. A 
pronounced Friedel-oscillation-like dependence of the coupling magnitude on the impurity position 
within the nanoribbon resulting from quantum size effects is found and investigated. In particular, 
the distance dependence of coupling is analysed. For semiconducting nanoribbons, this dependence 
is exponential-like, resembling the Bloembergen-Rowland interaction. In particular, for metallic 
nanoribbons, interesting behaviour is found for finite length systems, in which zero-energy states 
make an important contribution to the interaction. In such situation, the coupling decay with the 
distance can be then substantially slower. 
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I. INTRODUCTION 

Graphene, being an unique two-dimensional mate- 
rial 0, 0, opens new possibilities for contemporary 
physics and technology. One of the challenges is finding 
route to the spintronic applications of graphene and its 
derivatives 0, what motivates experimental and theo- 
retical studies of their magnetic properties. In particu- 
lar, interactions between magnetic moments introduced 
to graphene lattice focus considerable interest. One of the 
known mechanisms of interaction between magnetic mo- 
ments embedded in the host material is the Ruderman- 
Kittel-Kasuya-Yosida (RKKY) mechanism 0-0] of in- 
direct coupling mediated by the charge carriers. The 
form of RKKY interaction, originally derived for bulk, 
metallic, three-dimensional systems, is sensitive to the 
geometry and dimensionality of the underlying system 
and closely connected with the electronic structure of the 
host material. This opened room for studies of RKKY 
coupling in low-dimensional systems [cj-fl5l|. in particu- 
lar to mention thin films [Rj| - [l8| at the crossover be- 
tween the three- and two-dimensions, or for pseudo-one- 
dimensional systems [l9j . In particular, in metallic, two- 
dimensional systems, RKKY coupling has been predicted 
to decay with the distance according to J RKKY oc 1/r 2 . 
However, the peculiar dispersion relation for charge car- 
riers in monolayer graphene[2fjj. being a zero-gap semi- 
conductor, together with a bipartite nature of the under- 
lying graphene crystalline lattice with dominant nearest- 
neighbours electron hopping, causes the indirect RKKY 
coupling in graphene to differ from that found in 2D 
metals. The indirect RKKY exchange in graphene at- 
tracts considerable attention (2l| - |43j . It has been verified, 
that its asymptotic distance dependence is J RKKY cx 
1/r 3 for monolayer graphene (somehow like in a three- 
dimensional system). Moreover, in undoped graphene, 



the coupling between magnetic impurities in the same 
sublattice is always ferromagnetic, while for the impuri- 
ties in different sublattices, it is antifcrromagnetic. 

One of the promising building blocks for graphene 
spintronics are graphene nanoribbons. These nanostruc- 
tures offer highly tuneable electronic structure, vitally 
dependent on its geometry, mainly the edge form. In 
general, the form of the edge of graphene and similar 
structures is known from early experimental and theo- 
retical studies to shape electronic states in the vicinity 
of the Fermi level(e.g. [HI, El])- Besides zigzag-edged 
nanoribbons, offering expected spontaneous edge mag- 
netic polarization and thus concentrating the interest, 
also armchair-edged graphene nanoribbons (AGNRs) at- 
tract much attention (e.g. [46- 48]). The interest in prop- 
erties of AGNR is especially motivated by, for example, 
recent achievements in fabrication of narrow nanoribbons 
characterized by well-controlled width[49|, [5(|- On the 
other hand, the form of RKKY interaction in metal- 
lic nanoribbons with quadratic dispersion relation for 
the c harg e carriers has been also a subject of a current 
study [19] which emphasized the possibility of control over 
the coupling due to its position dependence. The last 
feature, resulting from the pronounced quantum size ef- 
fects, combined with an uniqueness of graphene, can be 
expected to promote interesting behaviour of RKKY in- 
teraction in graphene nanoribbons. Moreover, in gen- 
eral the effects of dimensional cross-over from quasi- 
one-dimensional nanoribbons to two-dimensional infinite 
monolayer appear interesting. This serves as motivation 
for computational studies of indirect coupling between 
magnetic impurities in armchair graphene nanoribbons 
presented in the paper. 



II. THEORY 
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The basis of the calculations is the choice of a finite 
graphene system. The study was performed for a finite 
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system, schematically depicted in Fig. [TJ The unit cell 
used in calculations contained N carbon atoms at the 
zigzag edge and the number of such zigzag atomic rows in 
armchair direction was equal to M, so in total the struc- 
ture contained M ■ N atoms. Empty and filled circles 
allow to distinguish between the sites belonging to two 
sublattices. The row of atoms extending along armchair 
direction is composed of subsequent nearest-neighbour 
pairs (with the nearest- neighbour distence equal to ao) 
and can be called a dimer line |5l|); the system contains 
thus N such dimer lines. The corresponding length of 
the nanostructure was L — 3aoAf/2, while its width was 
equal to W — V3ao (N — 1) /2. Periodic boundary con- 
ditions were used to connect the zigzag edges, so that for 
the zigzag atomic rows M + 1 = 1. The usage of such a 
geometry had two aims, first one being to simulate the 
behaviour of the infinite length AGNR, the second one 
being to investigate some important finite-size effects dis- 
cussed further. Only the nanoribbons for odd values of 
N were considered, since the interest was focused on the 
AGNRs possessing the symmetry axis. The symmetry 
axis of the AGNR was described by 5 = 0, while the 
armchair edges are shifted by <5 = ± (N — 1) /2 dimer 
lines from the center. 

The values of M used in the calculations of distance 
dependences of the coupling (Figs. 1315 1) were M = 402, 
while the rest of results, apart from those illustrated in 
Fig. 0Jb) was obtained using M = 202. It has been ver- 
ified that in absence of zero-energy states participating 
in indirect coupling (see further discussion), such a val- 
ues of M were large enough to guarantee the convergence 
of the results to the limit of infinite-length nanoribbon. 
Anticipating the further analysis, let us mention that the 
largest distance between the magnetic impurities consid- 
ered in this paper was t/ciq = 27, so that the padding 
around the impurity pair in the armchair direction was 
at least more than 10 times larger than that distance (see 
the results and discussion by Black-Schaffer[32J, demon- 
strating that the padding as twice as the distance be- 
tween the impurities is sufficient). 

The electronic structure of finite AGNRs was mod- 
eled using a tight-binding Hamiltonian for p z electrons, 
with nearest-neighbours hopping t, supplemented with 
an Anderson-Kondo impurity term [32[ : 

ft = -* ( c L c j> + c ]> c v) 




(1) 

where denotes summation over nearest neighbours 
in the lattice, while a = f, 1 is the electron spin. 

The coupling between on-site impurity spins S a and 
Sb (having the magnitude of S = 1 for simplicity and lo- 
cated at sites a, b) and the electron spins at the same sites 
is parametrized by the energy J, called further a contact 
potential. In all the further results, except those pre- 
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FIG. 1. Schematic view of the nanoribbon system used in the 
calculations. The solid box encloses the unit cell of width W 
and length L, with N carbon atoms at the zigzag edge. The 
dashed lines enclose a single row of atoms in zigzag direction; 
the unit cell is composed of M such rows. Periodic boundary 
conditions connect the zigzag edges, so that M + 1 = 1. The 
symmetry axis of the nanoribbon is marked with a dashed 
line. Open and filled circles denote the sites from two sublat- 
tices. 
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FIG. 2. Plot of the charge density associated with zero-energy 
states which are present in the electronic structure of metal- 
lic AGNR. Empty and filled circles depict the charge density 
for two distinct zero-energy states associated with two sub- 
lattices. The symmetry axis of the nanoribbon corresponds 
to S = 0, while at the armchair edges 8 = ± (N — 1) /2. 

sented in FigfJJa), the value of J/t = 0.1 was accepted. 
We note that the coupling is insensitive to the sign of 
J. Let us mention that for the sake of simplicity of nu- 
merical calculations, the Ising form of coupling between 
spins was assumed, and this assumption does not influ- 
ence the resulting value of indirect coupling integral be- 
tween impurity spins |4lj . The reduced Hamiltonian for 
impurity spin pair would be H a b = —J RKKY S^S^. In 
our convention the negative coupling values correspond 
to antiferromagnetic coupling, while the positive ones to 
ferromagnetic coupling. 

The Hamiltonian was diagonalized in single-particle 
approximation, with the help of LAPACK library [52j . 
yielding a set of electronic eigenenergies sorted in as- 
cending order. The ground state energy (at T = 0) for 
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the system depended on the orientation of the impurity 
spins and was expressed by a sum of M ■ N least eigen- 
values, E(S a ,Sb) — e i- This followed from the 
half-filling (charge neutrality) condition (one p z electron 
per each carbon atom) . Then, the indirect coupling con- 
stant J RKKY W as calculated from the expression [33 , l4l| : 

2S 2 jRKKY = E ( Sa = 1, S b = —1) — E (S a =l,S b = l). 

(2) 

Note that the approach used allows for non- 
perturbative determination of indirect exchange inte- 
grals. 

Let us clarify here that we refer to the resulting indi- 
rect coupling as RKKY interaction. In general, as shown 
further, the coupling may contain a highly significant 
contribution resulting from the mechanism which can be 
explained by means of first order perturbation calculus, 
in contrast to the original RKKY mechanism involving 
a second-order process. However, for brevity we call the 
total resulting coupling the RKKY coupling [4lj. 

We mention that due to lack of free zigzag edges in the 
assumed geometry of AGNRs, no additional zero energy 
states localized at zigzag edges are introduced this way 
to the electronic spectrum. 

The electronic structure of the AGNRs within nearest- 
neighbour tight-binding approximation has been exten- 
sively studied in the works of Wakabayashi et al. [53l — 
|{55j |. One of the results was that the AGNRs can be 
either metallic (M AGNRs), for widths N = 3n - 1 
(n = 1,2,...), or semiconducting (SC AGNRs), for 
N = 3n and N = 3n + 1. The SC AGNRs have an 
energy gap inversely proportional to the width. This be- 
haviour was fully confirmed in our calculations. Let us 
note that the two series of SC AGNRs, 3n and 3?i + 1, 
are inquivalent, so that two SC AGNRs from different se- 
ries can display some differences in the properties other 
than resulting merely from the different gap value. Such 
effects were for example predicted for excitonic spectra 
[EH HI HI- As N -> oo, the properties of AGNRs tend 
to that of an infinite graphene monolayer, being a zero- 
gap semiconductor. Let us emphasize that in our study 
we focus our interest only on the AGNRs with symmetry 
axis, therefore, the allowed values of N are only odd. 



A. Zero-energy states 

Let us focus on metallic AGNRs. In absence of the 
magnetic impurities, the electronic structure of metallic 
ANGRs is characterized by the existence of two zero- 
energy states (each one with an additional two-fold spin 
degeneracy) . The distribution of the charge density asso- 
ciated with both states is depicted in Fig.[3]for an AGNR 
with N = 11. First let us note that each of the states in- 
volves one sublattice only. Then, it follows that the wave- 
function amplitudes of both zero-energy states vanish for 
specific dimer lines, namely for that characterized by the 
distance from AGNR center equal to <5 = 3d, d = 0, 1, 



The number of such sites in our finite system is MN/3. 
In particular, the wavefunctions always vanish along the 
AGNR symmetry axis. The partial charges at the re- 
maining MN/3 lattice sites of each sublattice are dis- 
tributed uniformly and the probability of finding an elec- 
tron is then equal to uq = 3/ (MN) for each of this sites, 
as follows from the normalization condition. Especially, 
the wavefunction amplitudes take nonzero values at the 
AGNR edge and for the neighbouring dimer line. The 
form of the electron wavefunctions for AGNRs has been 
also discussed recently by Sasaki et aL[H[ (see the charge 
density for Dirac singularity zero-energy states sketched 
in the Fig. 1(b) in the Ref. l58l however, for a nanoribbon 
with even N). The wavefunction behaviour is a result 
of the interference of the electron wavefunctions incident 
and reflected from the armchair edge. The presence of 
similar interference effects has been also confirmed exper- 
imentally by means of STM/STS studies of the graphene 
sheets in the vicinity of the armchair edge [59, 60]. 



III. RESULTS 

In this section we present the results of numerical cal- 
culations performed for our system. 

One of the consequences of the peculiar Dirac-like dis- 
persion relation for charge carriers in graphene mono- 
layer is that the asymptotic distance dependence of the 
RKKY cou pling between magnetic impurities is 1/r 3 (see 
e -g- [HI HH mf), contrary to 1/r 2 characteristic of two- 
dimensional metallic systems [9hTT|. On the other hand, 
for metallic one-dimensional metallic systems the asymp- 
totic dependence on 1/r is expected (§. l9l. H2I llH] . Let us 
then analyze the distance dependence trends for the in- 
direct coupling between magnetic impurities in AGNRs. 
The distance dependences are presented, for representa- 
tive cases of three AGNRs with the widths of N = 7, 
9 and 11, in the Fig. [3] The separate plots are pre- 
sented for two semiconducting AGNRs and for a metallic 
AGNR; on the other hand the coupling is also plotted 
separately for the impurities in the same sublattice and 
in different sublattces. Let us note that in graphene, due 
to the bipartite nature of its lattice, the coupling be- 
tween impurities in the same sublattice is always ferro- 
magnetic, while for the impurities in different sublattices 
it is antiferromagnetic.!37] This rule is also obeyed in the 
case of AGNRs and no terms are present in the Hamil- 
tonian to cause its breaking. The coupling energies are 
depicted in semilogarithmic scale to illustrate better the 
dependence ranging over many orders of magnitude, thus 
the plots present the absolute values of the couplings. 
In all the plots, the dashed lines present the results of 
the perturbative calculation of the RKKY coupling en- 
ergy for armchair direction in monolayer graphene for 
J/t = 0.1 performed using the Eq. (19) and Eq. (26) in 
the Ref. [28| (for the same and for different sublattices, 
respectively). 

First let us observe that for semiconducting AGNR 
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FIG. 3. The distance dependence of an indirect RKKY coupling between magnetic impurities situated in different sublattices 
(left) and the same sublattice (right), for two classes of SC AGNRs [(a)-(d)] and M AGNR [(e), (f)]. Dashed lines denote the 
results of perturbative calculation after the Ref. l2g| . The plots are presented for various distances of the impurities from the 
AGNR center, covering the whole range between the center and the edge of the nanoribbon. 



with N = 7, the coupling tends to vanish with the dis- 
tance between impurities much faster than 1 /r 3 expected 
for graphene monolayer. This can be observed both for 
different sublattices (Fig. [UJa)) as well as for the same 
sublattice (Fig. GUb)). The dependence shows signs of 
exponential decay, according to J RKKY oc exp(— r/R). 
This kind of behaviour is observed for every distance of 
the impurity pair from the AGNR center, i.e. for every 
5; only some magnitude variations are present. Such an 
exponential decay is characteristic of the systems which 
indicate a gap in the electronic spectrum. For exam- 
ple, it has been mentioned in the Ref. [34] that the gap 
opening in graphene due to spin-orbit interactions should 
lead to such a behaviour. This kind of coupling is also 
known from the semiconductor physics as Bloembergen- 
Rowland coupling [6l| mediated by the excited electron- 



hole pairs. The characteristic decay distance R is energy 
gap-dependent and faster decay is predicted for narrower 
AGNRs with the larger energy gap. For the next width, 
N = 9 the AGNR is also semiconducting. Despite the 
small difference in the energy gap value, the distance be- 
haviour of the coupling is noticeably different. Namely, 
the decay with the distance appears significantly slower. 

The situation is very different for the metallic AGNR 
(studied case of N = 11). From the plots Fig. [3je),(f), 
it can be concluded that the coupling behaves in a man- 
ner characteristic of semiconducting systems only for the 
impurities distance from the AGNR symmetry axis 6 
such that S mod 3 = 0. For the remaining distances 
5, including especially the impurities at the AGNR edge, 
the distance dependence is much weaker. As it can be 
observed, for different sublattices, the AF interaction 
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tends to vanish quite slowly, and it can be numerically 
checked that it follows the J RKKY oc 1/r rule. Such a 
dependence is rather characteristic of metallic-like, one- 
dimensional systems. [§l-ficl [l3j . The distance dependence 
is even more interesting for the impurities at the same 
sublattice. Let us observe that the interaction energy, af- 
ter some initial decrease, is almost distance- independent. 
Moreover, the coupling magnitude for this case is greatly 
pronounced in comparison with the other cases. 

This peculiar behaviour of the indirect coupling can be 
understood owing to the analysis of the zero-energy states 
wavefunctions for AGNRs. Let us consider first the case 
of two magnetic impurities in different sublattices, situ- 
ated such that S mod 3^0. In such a situation, each 
one of the two zero-energy states is occupied by a sin- 
gle electron, with the spin orientation which minimizes 
the value of the term SkS^. present in the Hamiltonian 
Eq. (JXJ) . The first-order perturbation calculus correction 
to the energy of each state comes from the interaction 
with one impurity only and is equal to AEq = — ^uqS | J|, 
thus does not depend on the relative orientation of both 
impurity spins. Therefore, in such a situation zero-energy 
states give no first-order contribution to the indirect cou- 
pling. Let us also note that the zero-energy state is two- 
fold degenerate (not including spin degeneracy). The 
matrix elements of the Anderson-Kondo term taken be- 
tween two different zero-energy states vanish. As a con- 
sequence, the basis of the zero-energy states associated 
each with different sublattice constitutes a proper ba- 
sis for first-order perturbation calculus including the de- 
generacy. The second-order contribution emerges due to 
non-vanishing matrix elements of the Anderson-Kondo 
term taken between one of the zero energy states and 
other states. 

The situation is quite different when both magnetic 
impurities are situated in the same sublattice. If the 
impurity spins are aligned ferromagnetically, the zero- 
energy state belonging to this sublattice is occupied by 
a single electron with the spin direction which mini- 
mizes its energy. The remaining state (belonging to the 
other sublattice) is also occupied by a single electron, 
but it is unperturbed and its energy is unchanged. As 
a consequence, the first-order correction to the energy is 
AEq = — no'S'lJI. On the other hand, it the impurity 
spins are polarized antiferromagnetically, the first-order 
corrections to the energy of the state associated with 
their sublattice cancel each other. The other zero-energy 
state is still unperturbed, and therefore AEq F = 0. 
Taking in consideration the formula Eq. @, we ob- 
tain the indirect exchange integral associated with zero- 
energy states (in first-order perturbation calculus) equal 
to J RKKY = (AE AF - AE F ) /2S 2 , what leads to 
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RKKY 



= n \J\/2S. 



(3) 



Let us emphasize that this kind of contribution to the 
indirect coupling is absent when any state is filled with 
two opposite-spin electrons, since then the first order cor- 
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FIG. 4. (a) The dependence of the indirect coupling magni- 
tude on the value of contact potential J, for the impurities 
being 3rd and 4th neighbours, placed in the center or at the 
edge of a M AGNR with iV = 11. (b) The dependence of the 
indirect coupling magnitude on the length-to-width ratio for 
M AGNR with N = 11, for the impurities being 3rd and 4th 
neighbours, placed in the center or at the edge. 



rection to the energy of the state cancel out, while the sec- 
ond order corrections give rise to an ordinary RKKY in- 
direct coupling mechanism, with interaction proportional 
to J 2 . The first-order mechanism bears some resem- 
blance to double exchange and is somehow analogous to 
the ferromagnetic contribution to indirect coupling found 
in ultrasmall graphene nanoflakes [4l|. Due to the pro- 
portionality to | J\, it strongly dominates over the typical 
second-order mechanism for weaker contact potentials J. 
It is worth emphasizing that, despite the linear depen- 
dence on the contact potential, the coupling is insensitive 
to the sign of J. 

In order to illustrate this, let us study the dependence 
of coupling magnitudes on the interaction energy J be- 
tween the localized spins and charge carrier spins. The 
results of the calculations are depicted in double logarith- 
mic scale, in the Fig.@|a), for a selected metallic AGNR 
with N = 11. The magnetic impurities were placed ei- 
ther at the edge or on the symmetry axis and belonged 
to the same or different sublattices. For the impurities 
in different sublattices, it is visible that the coupling en- 
ergy J RKKY is proportional to J 2 . This means that the 
indirect interaction can be described using the second- 
order perturbation calculus (which is typically used to 
study indirect, charge-carrier mediated coupling in vari- 
ous systems). The situation is different when the impu- 
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FIG. 5. Distance dependence of indirect coupling magnitude for 4 selected metallic AGNRs (of widths from N = 5 to N = 23), 
for impurities in different sublattices (left) and in the same sublattice (right), for impurities placed in the center and at the 
edge of the nanoribbon. Dashed lines denote the results of perturbative calculation after the Ref. [28|. The plots are presented 
for various distances of the impurities from the AGNR center, covering the whole range between the center and the edge of the 
nanoribbon. 



rities belong to the same sublattice. Then, the quadratic 
dependence of J RKKY on J is conserved for the impu- 
rity pair in the AGNR center. In contrast, the interac- 
tion between edge impurities shows linear dependence of 
magnitude on \J\. This indicates that such an interaction 
contains a dominant contribution describable by means 
of first-order perturbation calculus, as it has been pointed 
out above with reference to the zero-energy states. The 
interaction magnitudes agree with those calculated from 
the Eq. J3]). In the studied range of J values, the cou- 
pling between the edge impurities is at least an order of 
magnitude stronger than for center impurities. 

It is of particular importance that the first-order con- 
tribution to the coupling is inversely proportional to the 
number of carbon atoms in the AGNR, thus it vanishes 
for infinite length of the nanoribbon and essentially con- 
stitutes a finite-size effect. This is a result of existence 
of just two zero-energy states for the whole system, for 
which n = 3/ (M • N). 

The importance of the finite length of AGNR for the 
mentioned ferromagnetic contribution to the coupling 
can be studied on the basis of Fig. HJb). This figure 
presents the indirect coupling magnitude as a function 
of L/W ratio for metallic AGNRs having the width of 
N = 11, for impurities in the same or in different sublat- 
tices, in AGNR center or at the edge. It is clear that for 
centrally placed impurities and as well for edge impurities 



in different sublattices, the coupling energy is indepen- 
dent on AGNR length for L/W ratios > 10 (which means 
achieving the limit of an infinite system size). In con- 
trast, for edge impurities in the same sublattice, the de- 
pendence of the type J RKKY cx 1/L is clearly observed. 
Note that even for L/W over 100, the coupling for edge 
impurities exceeds the one for center impurities by ap- 
proximately 2 orders of magnitude. Such a slow decrease 
establishes the importance of this contribution to indirect 
RKKY coupling. 

The evolution of distance dependence of indirect cou- 
pling when the metallic AGNR width is increased can 
be studied in Fig. [5] For impurities located in the cen- 
ter of the AGNR (Fig. EJa),(b)), the coupling decay is 
faster than 1/r 3 for all the illustrated AGNR widths, but 
the slope decreases gradually tending to reproduce the 
behaviour characteristic of an infinite graphene mono- 
layer. This kind of evolution can be expected for semi- 
conducting systems with a gap decreasing with the in- 
crease of the ribbon width. Recall that for the AGNR 
center, the zero-energy states are not involved in an in- 
direct coupling. Totally different situation is met for 
the impurities at the AGNR edge, for which the posi- 
tion zero-energy states can contribute. For edge impu- 
rities in different sublattices (Fig. |SJc)), the 1/r depen- 
dence is characterisatic of the thinnest AGNR case (what, 
as mentioned, resembled RKKY coupling behaviour for 
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FIG. 6. The coupling between (a) 3rd neighbour impurities 
(b) 4th neighbour impurities as a function of the AGNR width 
N. The coupling magnitudes are presented for the center 
and edge of the AGNR. Values for metallic AGNRs and two 
kinds of semiconducting AGNRs are depicted with different 
symbols. 



for J ft = 0.1 on the nanoribbon width (proportional to 
the number of atoms at the zigzag edge N). Fig. [5] shows 
the coupling between the impurities being the 3rd neigh- 
bours, situated either on the symmetry axis or at the 
AGNR edge. It is visible that a pronounced dependence 
of the coupling energy on the AGNR width occurs, with 
an oscillatory envelope. The oscillations tend to vanish 
for wider AGNRs. A period equal to 3 is visible in the 
changes vs. N, what corresponds with the width de- 
pendence of the electronic structure of AGNRs, as stud- 
ied using the tight-binding model by Wakabavashi (53h- 
[55j . The pronounced quantum size effects in RKKY 
coupling magnitude follow this rule. It can be also de- 
duced from the calculations for further neighbours that 
the oscillations envelope vanishes slower with the width 
for more distant impurities. Fig. [Sfb) depicts the case 
of the impurities being the 4th neighbours, at the edge 
or in the center of the AGNRs. The quantum size ef- 
fects for the impurity pair in the center resemble much 
ones observed for the previous case, with the same pe- 
riod of 3. The situation is much different for the edge 
impurities. There, a strong enhancement of the ferro- 
magnetic coupling is observed each time when N corre- 
sponds to the metallic AGNR. The origin of this coupling 
has been previously discussed. Let us observe that the 
energy of this coupling for metallic AGNRs varies ac- 
cording to J RKKY ex 1/N + const., in agreement with 
the described mechanism. The quantum size effects for 
AGNRs were already studied for example for the exci- 
tonic spectra [51], [56], [57j ■ 



one-dimensional metallic systems). It can be verified by 
plotting the dependence in double logarithmic scale (not 
shown) that the power decay law of coupling takes place 
for all studied widths of AGNR, however, the exponent 
varies from -1 toward -3 when N increases. Therefore, 
the distance behaviour of RKKY coupling evolves toward 
the expected behaviour for infinite graphene monolayer, 
for which J RKKY cx 1/r 3 is predicted. Moreover, the 
magnitude of indirect coupling is significantly reduced as 
a result of increase of the width. On the other hand, 
when the edge impurities are located in the same sub- 
lattice (the case illustrated in Fig. [5jd)), the influence 
of the AGNR width on the coupling magnitude is most 
pronounced. Under such conditions a first-order contri- 
bution to the coupling from a zero-energy state is present, 
but its value decreases inversely proportionally to the 
AGNR width (i.e. to the number of the lattice sites in 
the AGNR whilst the length is kept constant). For all 
the studied AGNR widths, the distance dependence of 
the coupling is very weak for larger distances. 

In order to analyse the width dependence on the cou- 
pling for the fixed distance between the impurities, let us 
focus on the case of the 3rd neighbours (situated in dif- 
ferent sublattices) and the 4th naighbours (the same sub- 
lattice). These cases were illustrated in the Fig. [13 which 
presents the dependence of the coupling energy J RKKY 



Finally, let us concentrate on the variations of the 
coupling magnitude for selected impurity pairs for var- 
ious pair distances S from the AGNR symmetry axis. 
Such calculation results are presented in Fig. [7J For the 
1st neighbours [Fig. [71(a)], the coupling behaves similarly 
for metallic and semiconducting AGNR, indicating pro- 
nounced oscillations while moving away from the symme- 
try axis towards the edge. In the vicinity of the edge the 
variations gain particularly large magnitude. This distri- 
bution of the RKKY coupling energies across the width 
of an AGNR is similar to Friedel oscillations, which are 
ubiquitous in numerous ph ysical p roperties of the sys- 
tems with boundaries (e.g. |62l |63|). For the 2nd neigh- 
bours [Fig. Efb)], the difference between metallic and 
semiconducting AGNR is clearly emphasized. For SC 
AGNR, the dependence of coupling magnitude resembles 
the Friedel oscillation-like picture from Fig.[7Ja). On the 
contrary, due to the existence of the mechanism enhanc- 
ing the ferromagnetic coupling with the help of the zero- 
energy state in M AGNR, the interaction is strongly pro- 
nounced for 5 ^ 3d (as explained before). Thus, the cou- 
pling energy is greatly pronounced except at that atomic 
rows for which the zero energy state charge density van- 
ishes. Let us note here that this enhancement would 
decrease in magnitude for longer AGNRs. 
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FIG. 7. The dependence of the indirect coupling on the dis- 
tance of the impurity pair from the AGNR center, plotted for 
the impurities being 1st neighbours (a) and 2nd neighbours 
(b), for two selected AGNRs - a SC one and a M one. 



IV. FINAL REMARKS 



gin, its distance dependence is predominantly of expo- 
nential kind, similar to Blocmbergen-Rowland coupling 
in gapped systems. This kind of behaviour appears also 
for metallic AGNRs when zero-energy states do not give 
contribution. 

On the contrary, the first-order mechanism becomes 
operative for the electronic states occupied only by a sin- 
gle electron, which can be formed in metallic nanorib- 
bons. Due to a peculiar form of the wavefunction for 
this states, the leading contribution to indirect coupling 
can be either slowly decaying (one-dimensional metallic- 
like in character) for different sublattices, or constant for 
the same sublattice. The mentioned effects are inevitably 
connected with a finite size of the system since their pres- 
ence is owing to the two zero-energy states only. How- 
ever, even for considerably large length to width ratio for 
AGNR, they could give a noticeable contribution to the 
coupling. It could be worth noticing that the involve- 
ment of zero energy states in charge carriers mediated 
coupling implies the large sensitivity of the interaction 
to the number of charge carriers in the system. This fea- 
ture in principle might open the door to control over the 
coupling sign and magnitude in AGNRs and may stimu- 
late further studies of such systems. On the other hand, 
also the influence of external fields on the system might 
be of interest. 

Finally, let us note that the presence of zero-energy de- 
generate states is characteristic of a wider class of carbon 
nanostructures (e.g. (64[) and thus such effects may not 
be limited to the system studied here. 



In the paper, the indirect (RKKY-like) coupling be- 
tween on-site magnetic impurities localized in an arm- 
chair graphene nanoribbon has been studied. The non- 
perturbational method, based on the total energy calcu- 
lation, allowed for capturing not only the usual RKKY 
mechanism attributed to second-order perturbation cal- 
culus, but also a first-order mechanism. It is worth 
emphasizing that the magnitude of the indirect RKKY 
coupling between magnetic impurity spins in AGNRs is 
strongly dependent on the location of the impurities in 
the nanoribbon, and the coupling energy exhibits pro- 
nounced oscillatory behaviour resembling the Friedel os- 
cillations. When the indirect coupling is of typical ori- 
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